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Abstract — We give analytical solutions to the titanium dioxide 
memristor with arbitary order of window functions, which 
assumes a nonlinear ionic drift model. As the achieved solution, 
the characteristic curve of state is demonstrated to be a useful tool 
for determining the operation point, waveform and saturation 
level. By using this characterizing tool, it is revealed that the same 
input signal can output completely different u—i orbital dynamics 
under different initial conditions, which is the uniqueness of 
memristors. The approach can be regarded as an analogy to using 
the characteristic curve for the BJT or MOS transisitors. Based 
on this model, we further propose a class of analytically solvable 
class of memristive systems that conform to Abel Differential 
Equations. The equations of state (EOS) of the titanium dioxide 
memristor based on both linear and nonlinear ionic drift models 
are typical integrable examples, which can be categorized into 
this Abel memristor class. This large family of Abel memristive 
systems offers a frame for obtaining and analyzing the solutions 
in the closed form, which facilitate their characterization at a 
more deterministic level. 

Index Terms — memristor, nonlinear dynamics, integrable sys- 
tem, device characterization. 

I. Introduction 

MATHEMATICAL tools have been inspiring and 
pedestalling the development in fundamental electron- 
ics, in parallel with experimental means, ever since the birth 
of this scientific branch. Memristors were first postulated 
theoretically by L. O. Chua in 1971 ingeniously on the 
observation of formal symmetries in the mathematical links 
among four elementary electronic quantities, and was proved 
on the ground of electromagnetics! 1]. The concept was re- 
vived by the realization of the titanium dioxide memristor 
in 2008|3|, which has been responded by massive studies. 
Consequently, the mathematical structure of passive electronic 
elements has been extended to a large family consisting of 
memristive, memcapacitive and meminductive systems|l2]||4l. 
And recently, the class of heterogeneous memristive systems 
has further broadened the concept from the prevailing homoge- 
neous systems|5|. Despite of the richness in the structures, the 
understanding of memristive systems is still far less accurate 
and systematic, as compared to the existing theories of BJT 
or MOS transistors. This is mainly caused by the lack of 
proficient tools in solving the governing nonlinear differential 
equations. 

In fact, since the emergence of the titanium dioxide mem- 
ristor, the dynamics of various systems have been mainly 
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Studied by numerical simulations. This is due to the rar- 
ity of integrable models (models that are solvable in the 
closed form) in memrisitve systems, unlike the cases of 
the various long-existing transistors, governed by algebraic 
equations. Numerical simulations have the merits of flexibility 
and versatility, but also suffer from the drawback that the 
physical relations of dominating factors to the dynamics can 
only be determined empirically or remain undetermined. In 
contrast, solutions in the closed form can offer the insight 
in such relations and render the behaviors reasonable and 
predictable. In this paper, we give the analytical solution to 
the titanium dioxide memristor with arbitary order of window 
functions, which takes into account the nonlinear ionic drift 
model. We demonstrate that the characteristic curve of state, 
as the achieved solution, is a useful tool in determining the 
operation point, waveform and saturation level, similar to 
the characteristic curve for BJT or MOS transisitors. At the 
meantime, a natural question is proposed whether there can 
be some classes of memristor models in general that are 
solvable or conditionally solvable in closed form. In Ref.||6l 
and Q, the authors heading for this question pioneered in 
proposing a class of memristive systems that conform to 
Bernoulli Differential Equations and argued that all current- or 
voltage-controlled memristors are classified therein. However, 
due to the implicit dependence of the memristance Ai{q{t)) 
(A4{cj){t))) on the current (voltage), this argument is only true 
for very special cases, e.g., the titanium dioxide memristor 
based on linear ionic drift model. We propose another class 
of (conditionally) solvable memristive systems that conform 
to (generalized) Abel Differential Equations. The equations 
of state of the titanium dioxide memristor based on both 
linear and nonlinear ionic drift models are typical integrable 
examples, which can be categorized into this Abel memristor 
class. 

II. Analytical Solutions to the Nonlinear Ionic 
Drift Model 

The titanium dioxide memristor, based on the mechanism 
of tunable doped region, was first implemented by B. Strukov, 
et al, and was formulated with a linear ionic drift modelOJ. 
However, the linear model fails in revealing the boundary 
effects, which should differ the rate of the change in the size 
of the doped region in the bulk from that at the boundary. 
To overcome this serious drawback, Y. N. Joglekar, et al, 
proposed a nonlinear ionic drift model in [81, which modifies 
the linear drift equation with a window function of defined 
orders p: 

F{x) = l-{2x-lf^' (1) 
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with p e M. It characterizes a milder drop in the rate of 
the dopant scale w at the boundaries with a smaller p, while 
a sharper drop with a larger p. The entire set of equations 
describing the dynamics therefore reads as 
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This is an analytically solvable model. By eliminating i{t) in 
Eq. (j2|, the equation of the state-variable w{t) can be written 

as 
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We name it as the Equation of State (EOS) of the memristor 
In fact, the EOS corresponding to different values of p can 
be categorized into a class of Abel Differential Equations |'9l. 
When p = 0, 1, Eq. ([3]) is an Abel equation of the second 
kind, which conforms to the equation 



W) + m ^ - f2{t)y{tf + h{t)y{t) + hit) (4) 

When p > 2, it is a generalized Abel equation in the form of 
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There have been solutions to vast forms of the Abel equations 
and well-established qualitative analysis tools. However, one 
can observe that the Abel equation in Eq. Q has a special 
form, which can be integrated directly by variable seperation. 
Making use of the special integral (for m G J\f) in the 
hypergeometric form 
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where 2^1 is the Gaussian hypergeometric function (also 
known as the (2, l)-ordered generalized hypergeometric 
function)! 12 1, the general solution to Eq. Q for any p E Af 
can be given in the following implicit form; 
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= fivRoN^{t) + C{y{0)) 

where denotes the integral of the input voltage signal 
(the flux) = J^v{T)dT, y{t) is a rescaled variable 

y{t) = 1 - 2w(t)/D e [-1, 1] and C{y{0)) is the integration 
constant determined by the initial condition y{0) (or w{0)). 
Equivalently, the solution can be also expressed in the form 
of series 
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by using the relation between the hypergeometric function and 
gamma functions, and noting the property T{z + 1) = zT{z): 
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In spite of the complexity in the general form, the analytical 
solutions to the system with window functions to the lowest 
orders can be given in a simpler form as follows, noting the 
Taylor series of logarithmic function and of inverse trigono- 
metric function arctan{x): 
p = Q : (with no window function) 
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where C{w{Q)) is the integration constant. The Abel dynamics 
of w(<) is the essential to the system, which links the voltage, 
current, flux and charge: 



q{t) ^ i{t) 
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III. Characteristic Curve for Determining 
Orbital Shapes 

Although the titanium dioxide memristor has been a focus 
of study in the area of memristive systems for long time, its 
dynamical behaviors still cannot been clearly characterized 
in a deterministic manner, compared to that for transistors. 
Numerical simulations, though flexible and versatile, are not 
proficient to reveal the inner relation between the parameters 
and the phase-planar dynamics. Critical questions, such as how 
the initial condition plays a role in the dynamics, which is 
crucial to a nonlinear system, and what determines exactly the 
deviation of the u — i orbital shape from a standard hyteretic 
loop, remain unexplained. The $ — w curve in the closed form, 
as given in Eq. ([8]), which we name as the Characteristic 
Curve of State, provides a useful tool for characterizing the 
dynamics of the memristor, which can determine the operation 
point, amplitude, waveform and saturation level, as using the 
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Fig. 1. Using the cliaracteristic curve "J" — to for cliaracterizing the 
Ti02 memristor of nonlinear ionic drift model, (a) and (b) demonstrate the 
determination of the operation point and the amplitude of the state-variable 
with the input signal and the initial conditions, under to(0) / D = 0.5 and 0.1, 
respectively (suppose p = 1). Q denotes the operation point, corresponding to 
the initial condition, (c) temporal variation in the memristance. (d) the u — i 
phase dynamics. The orbital shapes are determined by the waveform of w{t) 
in (a) and (b), including the saturations. 



characteristic curve in analyzing BJT or MOS transistors. The 
procedure is explained as follows: 

Suppose that the input voltage signal is v{t) ~ wq cos(a;i), 
which corresponds to a flux band with the width of 2<i>o = 
2vo/oj. For memristors, two factors are needed to determine 
the dynamical behavior of the nonlinear system: the input 
signal and the initial condition. We can see that the initial 
condition for the titanium oxide memristor actually functions 
as an operation point. The analytical solution in Eq. (j?]) 
indicates that different initial conditions C(z«(0)) are repre- 
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Fig. 2. Characteristic curve of transistors (BJT) as a tool for determining 
operation point, waveform and saturation level. The input cuiTcnt ii,{t) is 
transfered to the output current ic(t) and voltage Vce{t) with saturation at 
the peaks in the illustrated case. Q denotes the operation point along the load 
line. 



sented as a translation of the characteristic curve w{^), which 
determines the intersection of the characteristic curve and the 
input flux band. When the steeper part of the curve is within 
the interval [— $0i +*i'o]i the input signal $(t) is mapped onto 
the dopant scale w{t) with a larger amplitude; and vice versa. 
The amplitude of w{t) determines directly the change in the 
memristance A4{w{t)) (Fig. flic)) and hence the shape of 
the u — i orbital shape (Fig.^Td)). Similar to the case in 
transistors, the characteristic curve also reveals the saturation 
level in w{t), as it occurs, for example, in the upper peaks 
(w — )■ D) in Fig. ([T|a)). In fact, the saturation in w{t) will 
result in a deformation of the u — i orbital shape (see the 
triangular shaped orbit with a sharper negative resistance in 
Fig. [TJd) for w{0)/D = 0.5), as compared with a standard- 
shaped hysteretic loop (see Fig. [T|d) for w{0)/D = 0.1). It 
is of special interest that the same input signal can output 
completely different u — i orbital dynamics under different 
initial conditions, which is the uniqueness of memristors! 

The narrowing of the u — i orbital shape is usually ob- 
servable when raising the operation frequency IJI. However, 
the above analysis implies that this phenomenon can also be 
obtained by varying the initial condition. It is preceivable that 
a sufficiently large input signal will cause saturation on both 
upper and lower peaks (w — > and D) in w{t), which will 
produce a heavier deformation in the u — i orbit. The drastic 
difference in orbital shapes reminds us that to control the 
initial condition properly is critical in the correct operation 
of memristors, which is highly related to the reliability of 
memristive networks. 

In fact, this characterizing method can be well regarded as 
an analogy of using the characteristic curve in transistors for 
similar purposes (see Fig.|2]), in which the output signal ic{t) 
or Vce {t) encounters saturation when the operation point Q on 
the load line is close to the saturation region of the transistor 
Therefore, with the characteristic curve $ — w in the closed 
form, one is able to determine the main characteristics of 
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the dynamics of the memristor prior to amounts of numerical 
simulations in the same manner as for transistors. 

IV. Abel Dynamics in Memristive Systems 

As researchers have seen that analytical solutions can offer 
a deterministic characterizing tool for memristive systems, a 
natural question is raised whether there can be certain classes 
of analytically solvable memristive systems. This problem 
was first addressed by P. Georgious, et al. In Ref. |6| and 
an argument was made that any current-controlled mem- 
ristor(similar for the case of voltage-controlled memristor) 
belongs to the class of Bernoulli memristive systems. More 
specifically, the nonlinear first order ordinary differential equa- 
tion of the memristor on the current i{t): 
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is argued to conform to a Bemouli Differential Equation, in 
the general form of 
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where M.{q{t)) denotes the current-controlled memristance. 
As demonstrated in the referred papers|6||7|, this classification 
is correct for the titanium dioxide memristor based on linear 
ionic drift model as a current-controlled memristorp|. How- 
ever, this argument is not valid in general, due to the implicit 
dependence of the memristance M.{q{t)) (M{(j){t))) on the 
current (voltage). The titanium dioxide memristor model with 
nonlinear dopant drift, discussed in the last section, is a typical 
system beyond this category, despite of its current-controlled 
memristive nature. To see this point, it is noticeable that the 
relation between the charge q and the dopant scale w is 
modified by a hypergeometric factor for p > 0: 
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where B is an integration constant. Noting the relation in 
Eq. (9jl, q{w) is a series of odd-ordered terms of — l) 
and hence a bijective mapping between w and q (see Fig. 
|3]l. Therefore, the memristance AI{w) can be expressed by 
M[w{q)], which indicates that the system with window func- 
tions remains a current-controlled memristor. However, due to 
the relation of the term 
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in which dq{w)/dw is a function of w{t) by Eq. ( jTsj ), the 
right-handside of Eq. (13 1 has a dependence on the unsolved 
w{t). Henceforth, the u — i behavior of the system does not 
conform to a Bernoulli equation. 

Instead of Bernoulli dynamics, the EOS of the titanium 
dioxide memristor based on both linear and nonlinear ionic 




Fig. 3. The relation between the charge q and the dopant scale w. The 
p — > oo limit approaches the linear relation without window functions. 



drift models conform to the Abel Differential Equations, as 
demonstrated in the previous section. This categorization can 
be generalized to a large family of memristive systems- 
the Abel memristor class, which offers a large number of 
analytical solutions. Let us consider a memristive system with 
single state-variable. Suppose that v{t) is the input signal, the 
system equations are possible to be reduced to a differential 
equation on the state variable, by eliminating the current i{t). 
Assuming a general memristive system|4| 



v{t) = g{x, i,t) ■ i 
X = f{x,i,t) 
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can be reduced to a differential equation on the state-variable 
x{t) or more generally a bijective mapping of it y[x{t)], we 
define it as the general Equation of State (EOS) of the system: 



5(y(i),y(t),«(t),t) =0 
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This reduction is possible for most memristive systems, since 
the governing equations are algebraic equations of i{t). If the 
EOS has one of the following polynomial forms 

'^^^^^ - f3{t)yitf + h{t)y{tf + h{t)y{t) + m (19) 
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namely, the Abel Differential Equations of the first kind 
and of the second kind, we classify such systems into Abel 
memristive systems of the first kind and of the second kind, 
respectively. Actually, the Abel memristive systems of the 
second kind can be converted to the first kind under the trans- 
formation y{t) +i]{t) = 1/Y{t); inversely, an Abel system of 
the first kind can also be converted to an Abel system of the 
second kind under the transformation j/ — > y^^|9|. 

The Abel class is an important family in studying memris- 
tive systems for its formal commonly and solvability. Although 
the Abel equations are not integrable for arbitary coefficient 
functions, there have been solutions to a large family of 
specific forms, especially solutions to periodically driven 
systems. A standard procedure in solving an Abel equation 
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of the second kind is to make a reduction to the canonical 
form of n ■ riz — il, = R{z), with il^ being the derivative 
to the transformed variable z, to which a large number of 
solutions have been provided in Ref. ||9|- On the other hand, 
various tools for qualitative analysis have been established in 
analyzing the solutions and the stabilityllO|[ll |. 

The two classes of memristive systems can be further 
enlarged to the generalized Abel memristive systems with 
polynomials to higher orders 

dy{t) 



dt 



- fM(t)y(t)^' + ■■■ + f2{t)y{tY + h{t)y{t) + fo{t) 

(21) 



or its equivalent form 
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It is observable that when fk{t) 
1), Eq ^ 
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0, for k = 0,2, ...,(M- 
is again reduced to a Bernoulli Differential 
Equation|,9J (this is, however, an equation on the state-variable, 
instead of on the current or voltage). As supported by the 
achievements in the theory of ordinary differential equations, 
the generalized Abel class provides an important family of 
memristive systems in the polynomial form, under which a 
large number of analytical solutions are obtainable to the 
correponding models. 
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V. Conclusions 

We have given analytical solutions to the titanium dioxide 
memristor based on the nonlinear ionic drift model. The 
achieved characteristic curve of state as the solution has 
been demonstrated to be a useful tool in characterizing the 
memristor, similar to the characteristic curve for BJT or MOS 
transisitors. By using this characterizing tool, we have discov- 
ered that the same input signal can output completely different 
u — i orbital dynamics under different initial conditions, which 
is the uniqueness of memristors! Based on this model, we 
have further proposed the Abel class of analytically solvable 
memristive systems. The EOS of the titanium dioxide memris- 
tor based on both linear and nonlinear ionic drift models are 
typical integrable examples, able to be categorized into this 
Abel memristor class. This large family of Abel memristive 
systems has offered a frame for characterization of the systems 
at a more deterministic level via the solutions in the closed 
form. 
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